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We derive a dynamical mean-field theory for mixtures of interacting bosons and fermions on a lattice (BF- 
DMFT). The BF-DMFT is a comprehensive, thermodynamically consistent framework for the theoretical 
investigation of Bose-Fermi mixtures and is applicable for arbitrary values of the coupling parameters and 
temperatures. It becomes exact in the limit of high spatial dimensions d or coordination number Z of the 
lattice. In particular, the BF-DMFT treats normal and condensed bosons on equal footing and thus includes 
the effects caused by their dynamic coupling. Using the BF-DMFT we investigate two different interaction 
models of correlated lattice bosons and fermions, one where all particles are spinless (model I) and one 
where fermions carry a spin one-half (model II). In model I the local, repulsive interaction between bosons 
and fermions can give rise to an attractive effective interaction between the bosons. In model II it can also 
lead to an attraction between the fermions. 
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1 Introduction 

Ultracold atoms in optical lattices are not only fascinating many-body systems per se, they also help to 
better understand interaction models and correlation phenomena in condensed matter physics |[Tl|2l. In- 
deed, the high tunability of experimental parameters in these systems allows to explore the properties of 
correlated fermions and bosons in parameter regimes which are not accessible in condensed matter systems 
provided by nature Q. For example, precise measurements of the effects of strong disorder on interacting 
bosons have now become possible |[3|. Similarly, important insights can be expected from investigations of 
strongly disordered interacting fermions, for which detailed predictions exist |4, 5 , 6 , 7 , 8 , 9 , lOj. Ultracold 
atoms in optical lattices may also be employed to realize entirely new physical systems, e.g., mixtures with 
an arbitrary relative concentration of lattice fermions and bosons lfmfT2l[T3l . 

Quantum many-body problems in the thermodynamic limit can almost never be solved exactly and 
thus require approximate investigation methods. During the last 20 years the dynamical mean-field theory 
(DMFT) has proved to be a reliable approximation scheme for correlated lattice fermions (electrons) in 
dimension d = 3 lfT4l[T5l[T6l[r7ll . The bosonic counterpart of the DMFT — the B-DMFT — was developed 
only most recently lfT9ll . The B-DMFT is a thermodynamically consistent, conserving and non-perturbative 
theoretical framework for the investigation of correlated bosons on a lattice. It is applicable for all physical 
parameters and becomes exact in the limit of high spatial dimensions d or lattice coordination number Z. 
The B-DMFT has the important property that it treats normal and condensed bosons on equal footing and 
thus includes the effects caused by their dynamic coupling |fT9| . This is in contrast to the seminal static 
mean-field theory of Fisher et al. fZOl . and the equivalent Gutzwiller approximation 1211 . which focus on 
the condensed bosons and treat normal bosons as immobile. Explicit results of the B-DMFT were so far 
obtained for the bosonic Falicov-Kimball model ||T9l and the bosonic Hubbard model lfT9l l22l l23l 1241 . In 
this paper we extend the DMFT and B-DMFT and formulate a dynamical mean-field theory for mixtures 
of correlated lattice bosons and fermions. The theoretical framework derived in this way, referred to as 
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BF-DMFT in the following, is a generalization of both the fermionic DMFT and the B-DMFT. It provides 
a set of self-consistency equations which are valid for arbitrary coupling parameters and temperatures, 
and which describe the static and dynamic properties of a system of mutually interacting lattice bosons 
and fermions. These BF-DMFT equations are exact in the limit Z ^ oo and provide a comprehensive, 
non-perturbative approximation scheme for Bose-Fermi mixtures on lattices with finite Z. Similar to the 
B-DMFT the BF-DMFT derived here treats normal and (pair-)condensed particles on equal footing. 

Since the BF-DMFT includes many-body correlations and the dynamic coupling between normal and 
condensed particles in a systematic way it extends earlier theoretical investigations of Bose-Fermi mix- 
fiires formulated in the continuum Es] Eg] EH EH EH [H and on lattices 131] [32l [33l [34l [35l [36l [37l 
[38l |39l l40l l4n l42l l43l |44l |45l. A dynamic coupling between fermions in Bose-Fermi mixtures was in- 
cluded by Titvinidze, Snoek, and Hofstetter 1431 who treated fermions within DMFT and bosons within 
the Gutzwiller approximation. One of the most interesting questions in this context concerns the renor- 
malization of the interaction between bosons and between fermions due to the interaction between bosons 
and fermions. Employing the many-body BF-DMFT formulated in this paper we find that the effective 
interactions between bosons and between fermions can indeed be attractive, and can thus lead to phase 
instabilities towards superfluid fermions or phase-separated bosons. This confirms and extends the range 
of validity of earlier investigations performed within a T-matrix approximation or the static mean-field 
Bogoliubov approximation |E5lE6lE7l[3l][32lEll[321E3|39l|30l|4g|^ 

In section 2 we introduce two models of correlated, spinless lattice bosons interacting with correlated 
lattice fermions without spin (model I) and with spin (model II). In sections 3 and 4 we derive the BF- 
DMFT self-consistency equations for these models and discuss some of their properties, e.g., the emer- 
gence of effective attractions between bosons due to the interaction with fermions, and vice versa. 



2 Lattice models for mixtures of correlated bosons and fermions 

In the following we consider two different models for correlated lattice bosons (b) and fermions (/) with 
local interactions. In the first model (model I) all particles or atoms are spinless. Such a system can be 
obtained experimentally by preparing a mixture where ""^"K fermionic and ®^Rb bosonic atoms are in only 
one hyper-fine state ifTTl [121 [131 . This would correspond to the simplest Bose-Fermi mixture where all 
spin-dependent effects are absent. The Hamiltonian has the form 

^ = E ^Mb^ + E + Y E - 1) + E + E + f^/ E 

ij i ij i i 

where bi, fi (&|, f^) are annihilation (creation) operators on lattice site i for bosons and fermions, respec- 
tively, with number operators = blbi andnj ~ fjfi and hopping amplitudes andt^^ . The interaction 
between bosons is denoted by Ub, and that between bosons and a fermion by Ubf- In model I there is no 
local fermion-fermion interaction because of the exclusion principle B6l . For generality we include local 
potentials e^, e{, acting on bosons and fermions, respectively. They can either be applied externally such 
as by a harmonic trap, or act intrinsically as in the case of an additional random potential. 

In model II the bosons are still spinless, but the fermions now carry a spin one-half. Experimentally the 
^"K fermionic atoms should then be prepared in two hyperstates with equal population, whereas the '^^Rb 
bosonic atoms should remain in one hyperfine state. In such a situation superfluidity of condensed Cooper 
pairs also becomes possible. In this case the Hamiltonian takes the form 

= E4^1^.+E^«"^ + tE"^(-^-i) 

ij i 

ij(T i<7 <7 i(7 
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where a denotes the spin (or hyperfine state) of the fermions. In contrast to (1) model II contains also a 
local interaction U f between fermions with the opposite spins. 

We will now derive the self-consistent DMFT equations for model I in real space P9l . To this end we 
employ the cavity method lfT6l together with a scaling of the hopping amplitudes for fermions lfT4ll and 
bosons lfT9ll48l . 



3 Solution of model I within Dynamical Mean-Field Theory 

In the DMFT approach the local, site-diagonal Green functions are of particular interest. The local one- 
particle bosonic Green functions are defined as 



V {TTbl{T)bl{0))s, (r^6](r)6,(0))s. 



(3) 



where we use the Nambu notation to incorporate the off-diagonal, anomalous Green functions in the case of 
Bose-Einstein condensed bosons ifSOllSTI . and Tr represents the time ordering operator 152]. The Nambu 
spinor operators for bosons are defined by 

b,= f ) and hl^iblh). (4) 



To characterize the bosons completely we also have to consider the order parameter of the BEC phase, i.e., 
the condensate wave function, given by 

*,(r) = (b,(r))5.. (5) 
The square of the absolute value of *i(r) yields the number of bosons in the condensate at site i, i.e. 

nT"" = {h\[T))sAUr))s.=\\^^{r)\\ (6) 

The equilibrium condensate density is expected to be independent of time. This still allows for a time- 
dependence of the phase of $i (r) . The local, site diagonal one-particle fermionic Green function is defined 
by 

Gf(r)^-(T,/,(r)//(0)>5.. (7) 

The local functions defined above are determined by the local action Si which will be derived next. 

Using the path integral formalism ll53l the exact grand canonical partition function for the lattice Bose- 
Fermi mixture is obtained from 

Z = | D{h\ I D[f]e-^^'>'f\ (8) 
where the action S[b, f] is given by 

s[b, f] = jyri^J2 " ^'bMT) + f*{T){dr - ^i^)f^{T)] + Hi% /]| . (9) 

Here ^i, and are the chemical potentials of bosons and fermions, respectively, and (3 = l/ksT is the 
inverse temperature. In the cavity method lfT6l we select a single lattice site iq and split S[b, f] into three 
contributions: 

S[b, f] = S,, [6, /] + AS,„ [&, /] + S,^,, [b, /]. (10) 
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Here the first term. 



+ Y<W« W - 1) + U,MoiT)nliT) \ , (11) 



corresponds to the action for a separated site ig. The second term, 

AS,, [b, f] ^ f drY, {tl^bib, + i^^j^h, + ti^fy, + 4 /;/.„), 



(12) 



is due to hopping processes between site iq and all other sites, and the third term, Si^i, [b, /], represents 
the remaining parts of the action. In the next step we expand the exponential function with respect to 
ASig [b, f] and perform the functional integrals with respect to the complex variables 5; and Grassmann 
variables fi with i ^ iq. The resulting expression for the partition function contains infinitely many 
correlation functions with different numbers of operators. Within DMFT we need to keep only terms 
with a single bi, operator as well as pairs of bosonic or fermionic operators. These terms are then re- 
exponentiated according to the linked cluster theorem lfT6l[T9l . This approximation becomes exact in the 
limit of infinite dimensions d or coordination numbers Z, provided the hopping amplitudes for fermions, 
ifj 111411 and bosons, , |fT9l |48l are appropriately scaled with d orZ. In this limit the self-energies are 

diagonal in lattice indices, i.e., E^j'^(T) = T,^'-^ {T)Sij. As a result one obtains a sum of three local DMFT 
actions for the lattice site iq as 

S^o = Sl + Sl + S^. (13) 



Here 



= dThl{T){drCTs-{fi,-el)l)K,{T) + ^J^ dr j ^ drXMA^„(r-r')b,„(r') 

+ V r<M«(^) - 1) + /^^^ E «.o.b*„ (r)*,(r) (14) 
is the action for bosons coupled to a reservoir of normal and condensed lattice bosons lfT9l . 

Si = j'^ drfl (r) (dr -^ls+ eQ /.„ (r) + dr J^^ dr' (r) (r - r')/.o (r') (15) 
is the action for spinless fermions coupled to a reservoir of normal lattice bosons and 

St^^Utf drnl{T)ni{r) (16) 



"'0 

is the action describing the coupling between bosons and fermions at site iq. There is no coupling of 
fermions to condensed bosons since this coupling is described by disconnected correlation functions. The 
symbol 1 denotes a 2 x 2 unit matrix, and 0-3 is the Pauli matrices with ±1 on the diagonal. The factor 
1/2 is added because in the Nambu representation the number of degrees of freedom is doubled |50]. We 
will now discuss the physical meaning of the quantities A^^, A{^^ and Ki^j appearing in (fl4l i and ( fTsl l, and 
how they are determined. 

The quantities A^^ and A£ describe the resonant broadening of quantum-mechanical states on lattice 
site id. They may be interpreted as a hybridization of bosons and fermions, respectively, on site io with the 
surrounding bosonic or fermionic reservoirs and are given by 
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and 

A( (z^„) = - ^ t(,Gf^. (z^.„)4 . (18) 

Here Vn = 2n7r//3 and a;„ = (2n + l)Tr/P are even (bosons) and odd (fermions) Matsubara frequencies, 
respectively. The tilde denotes cavity Green functions which are determined for a lattice where the site 
iQ is removed. Nevertheless, using eqs. ( 117118b the hybridization functions can be expressed by local 
correlation functions only. Namely, since the self-energies are diagonal in the lattice indices and the cavity 
Green functions are related to the full Green functions by 

Gt.itiyn) = Gt^{^,y„) - Gi^{iv^)Gi,^{zv^)-'Gl^{zv^), (19) 



and 



one can derive p4| the corresponding local Dyson equations and express the hybridization functions by 
the self-energies and the local Green functions as 

Gliiv^r^ + = iv^a^ + (Mf, - <)1 - (21) 

and 

{il^n ) " ^ + S{„ (i CJ„ ) = iuJn + M/ - 4o - (*^" ) • (22) 

To close the set of the DMFT equations we employ the Dyson equation for the lattice Green functions, 
where the exact self-energies are replaced by the local ones, i.e.. 



and 



G{Aiuj„) = (iw„ + ^/ - lli{iujri))6ij - if- 



(24) 



The self-consistent solution is reached when the local Green functions Q and ^ are identical to the 
corresponding diagonal elements of the Green functions obtained from the Dyson eqs. (|23l ) and ( l24b . 

The scaled hopping amplitudes Ki^j appearing in eq. (fT4l i need to be determined such that in the non- 
interacting {Ub = 0) case the linear Gross-Pitaevski equation 

\ {drcr^ - (Mfc - ^^oir) + \ f dr' A^r - r')*,o(r') + ^ «.«*,(t) = (25) 

is satisfied for a given lattice ^55l . 

The real-space formulation presented here can be employed to study mixtures of bosons and fermions 
on a lattice in the presence of external potentials. It allows for the appearance of phases with diagonal 
long-range order, phase separation, and BEC. We note that when the external potentials are random, one 
also has to specify the average over the disorder |l5]|6l|2l. In the absence of external potentials (e^'^ ~ 0) 
there is usually no need for a formulation in real-space since the lattice Dyson equations can be written as a 
Hilbert transform over the momenta. However, even for e^'-^ = a real-space formulation is more suitable 
when one wants to describe phase separation or incommensurate long-range order, e.g., stripe formation. 
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3.1 Effective attraction between bosons mediated by their interaction with fermions 

Due to the absence of local fermion-fermion interactions the DMFT action for model 1 is bilinear in Grass- 
mann variables. Therefore they can be integrated out yielding an effective local action for the interacting 
bosons. The local partition function reads 



where the matrix elements of the operator M^^^ in Matsubara frequency space have the form 

L J nm 

iuJn - fJ-f + 4o + '^L (^") + ^"''o '^^^ ~ ^^^^ 

Apparently the local interaction between bosons and fermions Ub / gives rise to a complicated effective 
dynamics of the bosons. Indeed, even if the interaction between the bosons is neglected {Ub = 0) the 
effective hosonic action is not bilinear. We will now show that this corresponds to a nontrivial effective 
interaction between the bosons. To this end we expand the logarithm of the determinant as 

lnDet[M''] = Trln[A/''] = Trln[-(g^)"i + A/i] 

OO ^ 

= Trln[-(gO"'] - Y -TY[gfMlr, (27) 
^ — ' m 

m—l 

where we make use of the Weiss Green function 

Glo M = -— \ . . (28) 

By keeping only the first two terms one obtains the Random Phase Approximation (RPA), which takes into 
account bubble diagrams for fermions plus fermionic Hartree contributions. The effective bosonic action 
is then found as 

Si ^Sl + ^Y. si -(^)-^T.< ('^»)nfo i'^r.K (-^"), (29) 

where we introduced the fermionic polarization (bubble) function 

n{„ (l^n) = -4 E (^m + ^n). (30) 

The first-order term in [/f,y^ is the Hartree contribution and leads to a renormalization of the bosonic chemi- 
cal potential. The second-order term describes a frequency dependent (retarded) local interaction between 
bosons due to presence of the fermions. The static component of the polarization function is equal to 
the fermionic local density of states Nf^,^ at the Fermi level, i.e., n{^(;^„ = 0) = Nf^^{^). Therefore the 
effective static interaction between the bosons is given by U^^ ^ Ub — U^jN^^ (fi), which becomes attrac- 
tive for Ub < U^j:Nl^{fi). The bosonic subsystem then becomes unstable against phase separation or the 
formation of bosonic molecules. It is interesting to note that this effective attraction between the bosons 
originates from their interaction with the fermions. This is complementary to the well-known mediation 
of an effective attraction between fermions through the exchange of virtual bosons, e.g., phonons. This 
corroborates earlier results on the effective interaction between bosons Ii25. .26. .27. ,31. .32. ,28. ,37l |29l [39l 

Eiioiiiiiiiiiiiia. 
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3.2 Immobile fermions: Falicov-Kimball limit of model I 

In the case of immobile fermions, i.e., t{- = (atomic limit), model I reduces to a Falicov-Kimball like 
model 111911561 for Bose-Fermi mixtures. The number of fermions on each site is then conserved and can 
be used to characterize many-body states. The local partition function takes the form 

Z,„= e^^^f-^K f D\b\e-''o-^^^<^o'^'"^o^^\ (31) 

Similar expressions can be derived for bosonic correlation functions. Although the fermions are immobile 
they are still thermodynamically coupled to the bosons. This is expressed by the annealed (rather than 
quenched) average over different possible configurations of fermions on the lattice. In other words, due 
to the fluctuations around the microscopic equilibrium configurations the systems probes different con- 
figurations of fermions and, for a given temperature and other external parameters, realizes the optimal 
configuration. In particular, on a bipartite lattice and at low enough temperatures the fermions can be 
long-range ordered with a checker board structure. Therefore the t{j = (Falicov-Kimball) limit of model 
I differs from an Anderson disorder model, where fermions are fixed at random lattice sites which are 
independent of temperature or other parameters. 

For t{j = model I can have several different solutions. In particular, a sufficiently strong repulsion 
between bosons and fermions, Ubf, will cause a splitting of the bosonic band into two subbands. In this 
limit, and for a homogeneous system with boson density n'' satisfying nf' + nf G N, where is the density 
of fermions, a change of the interaction Ub will lead to a transition from a superfluid to a Mott insulator. 
Namely, for large values of Ub / the lattice sites occupied by fermions are inaccessible for bosons, thereby 
increasing the effective density of bosons. This transition is an analog of the metal-insulator transition for 
fermions at non-integer densities in the presence of binary-alloy disorder 1571 . Away from these special 
densities the ground state will be a superfluid. Since the condensation temperature Tb ec increases with 
the density of bosons we also expect Tbec to increase with increasing Ubf - Indeed, for [/;, = the model 
is equivalent to the bosonic Falicov-Kimball model with immobile, hard-core bosons discussed in Ref. ||T9l 
where it was shown that the transition temperature Tbec increases with increasing Ubf. 



4 Solution of model II within Dynamical Mean- Field Theory 

In the case of Bose-Fermi mixtures with spinful fermions described by model II, (2), condensation into 
a state with macroscopic quantum coherence is not restricted to bosons. Namely, fermions can undergo 
pair condensation provided ?7/ < 0, or if the effective boson-mediated attraction dominates over the bare 
repulsion Uf > (see below). To include the possibility of fermionic superfluidity (or superconductivity 
in the case of charged fermions) within the DMFT we employ the Nambu formalism also for fermions. To 
this end we define Nambu spinor fermion operators 

( J^T ) and f/ = (4,.fc) (32) 
and the local fermionic one-particle Green function 

\ [TrflAVfLmS, {Trflsij)hs{^))s, J 

For bosons all relevant Green functions and the BEC order parameter are the same as in model I. 

To derive the DMFT equations for model II in real space we again employ the cavity method and 
proceed as before. As a result we find the following local DMFT action at site 

S^,^S^,+Sl+S^l (34) 
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where 



Si - \[ drhl (r) {drCT; {fit - (r) + Ij^r J^^ dr^h^ (r) (r - r')b,„ (r') 

ni (r) K (r) - I) + ^ dr '^^o,K i^)^^ (^) (35) 
"'0 Jo , 



Ub 
' 2 



is the action for bosons coupled to a reservoir of normal and condensed lattice bosons |[T9l . 

Sfo = \[ ^^f^ (^) (^-1 - (a*/ - 4. >^ f^o (t) + dr'f * (r) a4 (r - r')f.o (r'), 

+ ^ rE<^M<.W (36) 
is the action for spinful fermions coupled to a fermionic reservoir, and 

St = f dr < (rhlAr) (37) 



is the action describing the coupling between bosons and fermions at site ig. The hybridization function 
matrices Aj^ and A{^ are related to the local self-energies through the local Dyson equations 

G,^„ (*i.„)-i + S?^ = z:.„a3 + (Mb - e,^Jl - A^„ (zz.„), (38) 



and 

GiiiuJn)-^ + = + (m/ - 4^3 - Af^(ic^„)- (39) 

The set of self-consistent DMFT equations for the model II is closed by the Dyson equations for the bosonic 
and fermionic lattice Green functions with local self-energies 

Gijitl^n) = [{iiyn(T3 + f^bl - ^^{^J^n))S^J - ^ l] , (40) 

and 



f ■ 



(41) 



These equations yield the general DMFT solution for Bose-Fermi mixture on a lattice. They include the 
possibility for bosonic and fermionic superfluidity, phase separation, and the appearance of phases with 
diagonal long-range order 

4. 1 Effective attraction between fermions mediated by the interaction with bosons 

In model II the presence of fermions with a local interaction can lead to similar effects as those discussed 
in connection with model I, namely to an effective attraction between bosons, and band splitting in the 
Fahcov-Kimball (atomic) limit where = 0. In particular, we will now discuss the situation where 
bosons influence the effective local interaction between fermions. For non-interacting bosons, Uh = 0, the 
local action is bilinear in the bosonic field which can therefore be integrated out analytically. We introduce 
a Bogoliubov shift of the complex bosonic field to remove the linear term in the bosonic part of the action 
as 

b{T)^b{T)+c{T), 5*(t)^&*(t)+c*(t). (42) 
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Here 6(r) and c(r) represent normal and condensed bosons, respectively. The local bosonic action then 
takes the form 



drbl (r) (a. - Mb + el + UtjnQ K (r) + I dr j dT% (r) (r - (r) 



JO 



(43) 



where the last term describes the Bose-Einstein condensate, and the chemical potential for bosons is renor- 
malized due to the interaction with the fermions by a Hartree-like term. The contribution describing the 
interaction between bosons and fermions is now expressed in terms of normal bosons b, i.e.. 



/3 



(44) 



In the next step we integrate out the normal bosons and thereby obtain the following partition function for 
fermions: 



Zl = / D 



■/]e-S4[/]-lnDot[M(J^ 



(45) 



Here the matrix elements of the operator Alf^^ in Matsubara frequency space have the form 



(dr + Ubfui + el + U,f J2 4j^rr' + A.^ (t - t') 



(46) 



Finally we introduce the Weiss Green function for local bosons 



(47) 



ii^n + M6 - Ubful - el - Al 
and perform a Taylor expansion to second order in Ubf- The effective fermionic action then takes the form 



of ^ of 



E^'o('^n)<(0) + ^ E <MTlli'^n)ni.'i-'^n). 



where the bosonic polarization function is given by 



(48) 



In the static limit this polarization function yields the local compressibility, n^^(t'„ = 0) = —dnl/dfit,, 
which must be positive for the system to be stable. The effective local interaction between fermions 
Uf — Ulfdn\ becomes negative (attractive) for Uf< U^fdn\ /dfib- Hence a repulsive in- 



eff 



Uy ' ~ u f — u f^j:U'ib'^^/ Ufn, uecomes negaiivc ^aiuacuve; loi u f ^ '^bf'^"'io/ 
teraction between fermions and bosons can lead to an effective attraction between fermions. The fermionic 
subsystem then becomes unstable with respect to a condensation of Cooper pairs as is well-known from 
the theory of superconductivity ll58l . 
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5 Summary 

We derived a comprehensive, thermodynamically consistent theoretical framework for the investigation of 
mixtures of correlated lattice bosons and fermions — the BF-DMFT. In analogy to its purely fermionic and 
bosonic counterparts the BF-DMFT becomes exact in the limit of high spatial dimensions d or coordination 
number Z. It may be employed to calculate the phase diagram and thermodynamics of mixtures of inter- 
acting lattice bosons and fermions in the entire range of microscopic parameters. As in the B-DMFT the 
BF-DMFT employs a different scaling of the hopping amplitude of bosons with Z, depending on whether 
the bosons are in the normal or the condensed phase. Thereby normal and condensed bosons are treated 
on equal footing. Namely, in the BF-DMFT the normal bosons retain their dynamics even in the limit 
Z ^ oo, such that the effects of the dynamic coupling between normal bosons and the condensate are 
fully included. 

Using the self-consistency equations of the BF-DMFT we solved two different interaction models of 
correlated bosons and fermions: model 1 where all particles are spinless, and model 11 where fermions 
carry spin one-half. In the latter model a local interaction between fermions is present which leads to 
dynamical effects even in the limit of large Z. In model 1 we showed that the local interaction between 
bosons and fermions can give rise to an effective attraction between bosons. The bosonic subsystem is 
therefore unstable against phase separation or the formation of bosonic molecules. In the case of immobile 
fermions model I reduces to a Falicov-Kimball model for Bose-Fermi mixtures which can be solved in 
closed form. Model II has several different solutions depending on the interaction parameters and the 
density of bosons. In particular, the interaction between bosons and fermions can mediate an effective 
attractive local interaction not only between bosons but also between fermions. The latter leads to an 
instability with respect to Cooper pair formation as in superconductivity. 

In conclusion, interaction effects in mixtures of bosons and fermions on a lattice can lead to a multitude 
of fascinating correlation phenomena, including the emergence of an effective attraction between bosons 
mediated by fermions and vice versa. Experiments with cold atoms in optical lattices will be able to test 
these predictions and lead to deep insights into the overall physical behavior of many-body systems with 
particles obeying different statistics. 
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